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Chiral phase transition in thermal QCD is studied by using the Dyson- 
Schwinger (DS) equation in the real time hard thermal loop approximation. 
Our results on the critical temperature and the critical coupling are signifi¬ 
cantly different from those in the preceding analyses in the ladder DS equation, 
showing the importance of properly taking into account the essential thermal 
effects, namely the Landau damping and the unstable nature of thermal quasi¬ 
particles. 


Although lots of efforts have been made to undestand the temperature/density-dependent 
phase transition in thermal QCD, we still have not really understood yet even the relation 
between the chiral and the deconhnement phase transitions. The analysis with the use of Dyson- 
Schwinger (DS) equation is a powerful tool to investigate with the analytic procedure the phase 
structure of gauge theories. At hnite temperature/density, however, we should seriously ask 
whether the preceding results [1] of the DS equation analysis on the chiral phase transition could 
be relied on being the real consequences of thermal gauge theories. This is because in previous 
analyses the lessons from vacuum theories have been applied without close examination, having 
missed the essential thermal effects, thus causing the neglection of large contributions to the 
DS equation unless otherwise existed. 

We here make a re-analysis by studying the hard-thermal-loop (HTL) resummed DS equa¬ 
tion in the real time formalism, giving a new understanding on the phase structure and the 
mechanism of phase transition in thermal gauge theories. Main interest of the present investiga¬ 
tion lies in clarifying what are the essential temperature effects that govern the phase transition 
and also in Ending how we can closely take these effects into the ’’kernel” of the DS equation. 
Essential procedures of our analysis can be summerised as follows; 

i) We use the real time closed-time-path formalism and study the physical mass function E/j 
itself, not the Sn, of the retarded fermion propagator. Sr is the mass function of ’’unstable” 
quasi-particle in thermal field theories, thus having non-trivial imaginary parts as well as non¬ 
trivial wave-function renormalization constants. 

ii) To closely estimating the dominant temperature effects we focus on studying the DS 
equation being exact up to HTL approximation: both the gauge boson propagator and the 
vertex functions are determined within the HTL approximation, guaranteeing the gauge invari¬ 
ance of the result at least in the perturbative analysis. With the HTL vertex functions [2] we 
can derive the desired HTL resummed DS equation. 

The second point listed above is better to be taken step by step into the actual analysis 
of the DS equation. In the present paper we present the result of our first step investiga¬ 
tion; focussing on what happens when taking into account exactly the HTL resummed gauge 


boson propagators. We show explicitly the results of preceding analyses can not be the real 
consequences of thermal gauge held theories. 

Decomposing the fermion mass function S/j as 

EniP) = (1 - A{P))pa^ - B{Ph^ + C{P) , (1) 

with A{P), B{P) and C{P) being the three independent scalar invariants to be determined, 
the DS equation we here investigate in the point vertex approximation with the exact HTL 
resummed gauge boson propagator becomes as (P^ = (po, p)), 
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where the cutoff scale A is introduced, regularizing the integral. *G'Y{K) in Eqs.(2-4) is the 
HTL resummed retarded gauge boson propagator [3], 
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where being the HTL resummed transverse/longitudinal mass function [4], and 

j^p(T _ gpa _ j^pa _ jgpcr^ j^pcr _ —YcpR'^/K'^^ DP"^ = RPR"^/K‘^, bciug the projection tensors. 
In the above, K = {k, /cok), k = and k = k//c. 



The DS equation, Eqs.(2-4), is still quite tough to be attacked, forcing us further approx¬ 
imations for the analysis to be effectively carried out. However, the approximation made use 
of must be consistent with the HTL approximation, without missing the important thermal 
effects out of the kernel of the DS equation. 

Here it is worth noticing that the instantaneous exchange (IE) approximation frequently 
used in the preceding analyses [1] is not compatible with the HTL approximation in the strict 
sense. In the exact lE-limit the HTL resummed transverse mass function, vanishes 

and the transverse (magnetic) mode becomes totally massless. Namely the IE approximation 
discards the important thermal effect coming from the Landau damping, thus dismissing the 
dynamical screening of the magnetic mode, causing the famous quadratic divergence of the 
Rutherford scattering cross section. This can be clearly seen by taking the lE-limit of the 
DS equation, Eqs.(2-4), and neglecting Im[A{P)] and Im[C{P)], then obtaining the following 
equation for Im[B{P)] {E = ^(p — k)^). 
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showing Im[B{P)] to be quadratically divergent. The reason why in the previous analyses this 
divergence problem did not appear is that the imaginary part of is completely neglected 
there from the beginning, namely the equation for Im[T,fi\ is totally discarded. 

Taking the above into account the approximation we further make use of is the IE ap¬ 
proximation to the longitudinal gauge boson propagator, by keeping the exact HTL resummed 
transverse propagator. In the lE-limit the HTL resummed longitudinal mass function, 
has a nonvanishing thermal Debye screening mass ~ thus even in the IE limit the 

longitudinal mode can take into account the essential thermal effect. In the present analysis 
the gauge is hxed to the Landau gauge (.^ = 0). 

Now we solve the DS equations Eqs.(2-4) with the IE approximation to the longitudinal 
mode, using the iterating procedure with suitable trial functions for the solution. This method 
is useful so long as the convergence of the iteration is guaranteed. 

At each iteration, the three-fold integration is performed over k, ko and = cos 6, 6 = p - h. 
The integration kernel of the present DS equation shows a little bit singular behavior, and the 
numerical integration of such a singular integrand needs careful integration prescription, which 
is properly managed. As a result we performed the 1000 times iterations, obtaining fairly stable 
solutions. 

Result of the present analysis shows that the wave function renormalization constants receive 
10-20 percent corrections, indicating the necessity of gauge-parameter dependent analysis. The 
generated size of the imaginary part is nearly the same as that of the real part, indicating the 
existence of the non-trivial imaginary parts, which is as expected though completely neglected 
in the previous analyses. 

Let us now show the behavior of the mass function at some hxed value of p as a function 
of the parameters a = g’^Cp/in and T. The T dependence of the mass function Re[C{P)] 
for various hxed values of a, and the a dependence of Re[C{P)] for various hxed values of T 
are shown in Figs.la and lb, both with po = 0 and p = O.IA. The errors resulting from the 
huctuations are smaller than the size of the symbol used for each sample point in Figs.la and 
lb. From these hgures we can see the two facts; i) The chiral phase transition is of second 
order, since a fermion mass is generated at a critical value of the temperature T or at the 
critical coupling constant a without any discontinuity, and ii) the critical temperature at 
hxed value of a is signihcantly higher than the previous results [1], namely the restoration of 
chiral symmetry occurs at higher temperature than previously expected. Our result shows that 



as T becomes smaller the critical coupling constant etc also becomes smaller, being consistent 
with the zero temperature result, where the critical coupling constant = tt/S. The second 
fact clearly shows that in the previous analyses the important temperature effects are neglected 
due to the inappropriate approximations. 

Present result shows the correctness of our research-strategy, thus the importance of the 
full HTL resummed DS equation analysis of the chiral phase transition at hnite tempera¬ 
ture/density, which is now under investigation. 
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